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The spatial patterns in species richness and abundance were investigated here at a 
series of reservoirs at different elevations, to establish which factors determine spe¬ 
cies distributions along this topographic gradient. Larvae of 18 spp. were sampled in 
small reservoirs across a 1250 m elevational gradient at one latitude. Most spp. oc¬ 
curred throughout all elevations indicating that this subtropical odon. assemblage as 
a whole is remarkably tolerant of elevational changes. Although Anisoptera larval 
species richness and abundance increased significantly with increasing elevation, there 
was no change in Zygoptera species richness, while Zygoptera abundance decreased 
significantly. Species-site-variable triplots for Anisoptera and Zygoptera larvae indi¬ 
cated that no measured site variable on an individual basis clearly accounted for lar¬ 
val species assemblage distribution patterns. Nevertheless, canonical axes and their 
respective intra-set correlation coefficients showed that some measured site variables 
e.g. floating/submerged vegetation, turbidity, pH, water temperature (resulting from 
sunny or shade habitat conditions), marginal grasses, water depth as well as elevation 
to some extent, explained the main variation in species assemblage composition/dis¬ 
tribution in a broadly similar manner for both suborders. Generally, the reservoirs 
recruited spp. from the regional pool, irrespective of the elevation of the pool. These 
spp. were all geographically widespread spp. that took advantage of the presence of 
these man-made reservoirs, and included only one national endemic. Although these 
artificial water bodies are not increasing the ‘extent of species occurrence’, they play 
a major role in increasing ‘area of occupancy’. Furthermore, these spp. are remark¬ 
ably vagile, habitat-tolerant, as well as elevationally-widespread. A reasonable expla¬ 
nation is that this assemblage is the historical survivor over many millennia of oscillat¬ 
ing wet/dry periods and natural selection. The present-day spp. are those that readily 
recolonise pools after drought has been broken, and are pioneering residents of new 
water bodies over a wide elevational range. 
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INTRODUCTION 

Reservoirs, being artificial, provide an increase in the area of occupancy for 
aquatic invertebrates. This increases regional invertebrate biomass which inevi¬ 
tably leads to increases in water bird and fish populations (WARD, 1984; Mc- 
CAFFERTY, 1984). Aquatic environments often contain large assemblages of 
co-existing resident Odonata species and more transient species (CORBET, 1999; 
STERNBERG, 1994). In small, still-water habitats, they are found mostly in the 
littoral zone (JOHNSON & CROWLEY, 1980; CROWLEY & JOHNSON, 1982), 
and play an important role in pond food webs as predators of other invertebrates 
(BENKE, 1976, 1978; THORP & COTHRAN, 1984; SAMWAYS et al., 1996) 
and even some vertebrates (CALDWELL et al., 1980; TRAVIS, et al., 1985), and 
in turn, are prey for fish (BENKE, 1978; THORP&BERGEY, 1981; JOHNSON 
et al., 1984; MORIN, 1984a, 1984b; PIERCE et al., 1985). 

Within the larger aquatic habitat, microhabitat preferences differ from one 
species to another, and even between larval instars of a species (CROWLEY & 
JOHNSON, 1982). Larvae may be cryptic or nocturnal but still suffer high natural 
mortality, so that the numbers of larvae in a given cohort decrease substantially 
as development proceeds. Also, as dragonfly larvae are predators, their numbers 
may be low (BENKE & BENKE, 1975; BENKE, 1976; MACAN, 1977). There 
are relatively little quantitative data on their population dynamics. Some studies 
describe patterns of larval abundance and survivorship for a single species (COR¬ 
BET, 1957; LAWTON, 1970;UBUKATA, 1981), and others provide information 
on a number of co-existing species in still-water habitats (BENKE & BENKE, 
1975; CARCHINI et al., 1992; WISSINGER, 1988). 

In South Africa, taxonomic difficulties have been a significant challenge in drag¬ 
onfly larval studies. The larvae of many species are unknown and virtually all 
early-instar larvae are unidentifiable. This means that it is difficult to determine 
the dragonfly assemblage in a habitat, let alone their abundance. 

Any changes in landscape will affect both the aquatic larva and the aerial adult. 
Little is known about which aspects of artificial water bodies encourage particu¬ 
lar species of African larval odonates either over space or time. The east coast 
of South Africa is interesting in this regard, because as well as having a major 
elevational gradient (0-3000 m a.s.l.) at one latitude along a 200 km transect, it is 
also subject to major changes in water on a periodic basis, as a result of El Nino 
Southern Oscillation (ENSO) events. The aim here was therefore to determine 
the rate of species turnover with elevation, and to establish which environmental 
factors might be responsible for any assemblage change. It was postulated that 
wide elevational and habitat tolerance was the end-product of considerable fluc¬ 
tuations in weather and climate. 
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METHODS 


Five reservoirs, across a 1250 m elevational gradient (100-1350 m a.s.1. in KwaZulu-Natal, South 
Africa, were selected (Fig. 1, Tab. I). They were: SB = Stainbank Nature Reserve (100 m a.s.1.), KL = 
Krantzkloof Nature Reserve (450 m a.s.1), BG = Pietermaritzburg Botanical Gardens (790 m a.s.1.), 
CE = Cedara (1050 m a.s.1.) and GH = Goodhope Estate (1350 m a.s.1.). 

Each reservoir was stratified into six sub-sites (SS) measuring 20 m length (along the pond edge) 
by 2 m width (1 m on land and 1 m into water). Data were collected on 42 sampling occasions from 
January 2001-December 2002. Abundance of larvae was estimated at each SS between 1 lhOO and 
13h00 using a dip-net (41 cm diameter x 1 mm mesh sieve). Two dips per SS (12 dips/site) were done 
within 20 min. Each dip was followed by vigorously shoving the net back and forth in the water and 
among water weeds, along rushes and beside banks. Larvae were identified using a hand lens and re¬ 
leased back into water. Where individuals could not be identified, they were retained for subsequent 
rearing and identification. Larvae of most species were reared in 2-3 cm depth of dechlorinated water. 
Larvae were then transferred from the collecting phials to small rearing cages made out of plastic net¬ 
ting rolled into cylinders of about 6 cm width and 15 cm height and sealed at the top with an inverted 
plastic petridish. To avoid cannibalism and to ensure correct association of larvae and adults, only one 
larva was placed in each cage. The water inside the cage was artificially aerated. Food was provided 
principally as mosquito larvae and young tadpoles netted from a nearby pool or rain barrel. 

Larvae were fed daily and uneaten or rotten food removed from the cage. After several hours, and 
usually during the night or early morning, the larva would leave the water and climb up the side of 


the cage. It usually stopped 
at the highest point where 
it emerged. The remaining 
exuviae was collected and 
retained in a phial contain¬ 
ing 75% ethyl alcohol. Af¬ 
ter a few days, the adult was 
identified relative to the ex¬ 
uviae. 

ENVIRONMENTAL 
VARIABLES. - Dis¬ 
solved oxygen and temper¬ 
ature influence the develop¬ 
ment of dragonfly larvae 
(CORBET, 1999), making 
shallow water particularly 
important for them. Many 
Odonata also prefer water 
of a certain pH (WEIR, 
1974; OSBORN, 1995). 
High turbidity may reduce 
visual detection of prey, 
and current rate may influ¬ 
ence oxygen concentration. 



Vegetation serves important Fig. \ Geographical location of study sites (KwaZulu, Southern Af- 
functions for larval Odona- nca): (A) Kenneth Stainbank Nature Reserve [SB] (100 m a.s.1.); - 
ta (BUCHWALD, 1992), (It) Krantzkloof Nature Reserve [KL] (450 m a.s.1.); — (C) Botanical 
and along with substrate, Gardens [BG] (790 m a.s.1.); - (D) Cedara [CE] (1050 m a.s.1.); - (E) 
provides refuge from preda- Mon di Goodhope Estate [GH] (1350 m a.s.1.). 
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Table I 

The five elevational sites used in this study 


Grid 

Site name/code reference Land use 

and Elevation 


Kenneth Stainbank 

29°50’S; 

Nature reserve 

Nature eserve (SB) 

30°55’E; 


(Low elevation) 

100 m a.s.1. 


Krantzkloof Nature 

29°46’S, 

Nature reserve 

Reserve (KL) 

30°5’E; 


(Mid-low elevation) 

450 m a.s.1 


Botanical Gardens, 

29°35’S; 

Botanical gardens 

Pietermaritzburg 

30°25’E; 


(BG) 

790 m a.s.1. 


(Mid-elevation) 
Cedara (CE) 

29°61’S; 

Low-intensity agricultural area 

(Mid-high elevation) 

29°06E; 


Mondi Goodhope 

1050 m a.s.1. 
29°40’S; 

Extensive grassland corridors (natural state) among 

Estate (GH) 

29°58’E; 

pine stands 

(High elevation) 

1350 m a.s.1. 



tors. The thermal requirements of the species results in different responses to sunlight and shade. 

Against this background, two one-monthly measurements of variables were made in sunny, warm 
conditions. Water depth was measured using a V.C 1456 meter rule. Acidity (pH) was measured using 
a Jenway 3405 electrochemistry analyser. Turbidity was estimated at midday as the results became less 
reliable near dawn or dusk because of reduced surface illumination. Percentage estimation was done 
visually, and ranged between 0% (totally transparent) to 100% (totally muddy conditions) using a 
Secchi disc. The degree of shading of SSs was estimated by mean percentage shade cover at midday 
for each SS. Water and ambient SS temperatures were measured using the Delta Trak hand thermom¬ 
eter at midday, two cm below the water surface. Ambient temperature records were cross-checked 
with the weather bureau. Rainfall and ambient and water temperature (At/Wt) data for Stainbank 
(SB) and Krantzkloof (KL) Nature Reserves were compared with that collated by the Durban Air¬ 
port weather station. Similar meteorological data for the Botanical Gardens (BG), Cedara (CE) and 
Goodhope Estate (GH), were compared with that collated by the weather bureau at Cedara Agricul¬ 
tural college. Sub-sites also classified on the amount of open water and whether aquatic vegetation 
covered the bottom of the water body. Vegetation was classified into: Submerged (Sveg), emergent/ 
floating (EFveg) and marginal (Mgra) vegetation and was estimated using percentages of sub-sites 
they covered. Other environmental variables measured were percentage of marsh (Marsh), Shade (% 
Sh) and elevation (Alt) (m a.s.1.) 

DATA ANALYSES. — Analyses of species spatial variability was done using linear regression 
and non-parametric methods. Spearman s correlation coefficients were used to measure the associa¬ 
tion between SS variables and species abundance and richness. These correlations were calculated us¬ 
ing SPSS version 6.1. GOTELLI & COLWELL (2001) have suggested that for patchy distributions, 
individual based rarefaction curves inevitably overestimate the number of species that would have 
been found with less effort. They recommend species accumulation curves that take into account be- 
tween-sample heterogeneity (COLWELL & CODDINGTON, 1994; UGLAND et al., 2003). Rela¬ 
tionships between abundance and the number of species having that abundance was shown graphi- 
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cally using plots of abundance of each species in the samples in order of its rank from the most to 
the least abundant for each site. This method presented a graphical and more explicit way of deter¬ 
mining species abundance patterns. 

Ordination describes a set of techniques in which sub-sites are arranged in relation to one or more 
coordinate axes such that their relative positions to the axes and to each other provides maximum 
information about their ecological similarities. When sub-sites that are most similar or dissimilar are 
identified based on coordinate positions, underlying biotic and abiotic factors that might be responsi¬ 
ble for the patterns are determined. The computer software CANOCO (TER BRAAK, 1986) which 
combines into one algorithm correspondence analysis (CA) on species data with weighted multiple 
regressions on environmental data was used. This technique relates species composition to known 
variations in the environment. CANOCO leads to an ordination diagram in which points represent 
species and sites, and vectors (arrows) represent environmental variables. Such a diagram shows the 
patterns of variation in species composition that can be explained best by the environmental vari¬ 
ables and also visualises approximate ‘centres’ of the species distributions along each of the environ¬ 
mental variables (TER BRAAK & LOOM AN, 1995). All abundance and environmental variable 
data were log-transformed to maintain normality and to satisfy the requirements of ANOVA and 
multivariate analyses. 


RESULTS 

OVERALL FAMILY/SPECIES TRENDS AT THE DIFFERENT ELEVATIONS 

A total of 18 species were sampled (Tab. II) alongside dominant subsite veg¬ 
etation made up of mostly indigenous plants. The elevational tolerance of larval 
odonates shows that most species, except Lestes tridens , occurred across a wide 
elevational band. An estimate of total species richness across sites using sub¬ 
site-groupings showed that both Anisoptera and Zygoptera had broadly similar 
trends (Fig. 2a, b). The asymptote was not reached for either suborder, despite 
much sampling, suggesting other species were still present but not sampled. On 
the whole, there was much evenness between abundance of species. At 100 m a.s.1. 
(SB), Anax imperator , followed by Crocothemis erythraea larvae were the most 
dominant Anisoptera species at 35% and 33% dominance levels respectively. The 
remaining elevation sites had high evenness and were similar in patterns of larval 
species dominance, ranging from 19% at 1350m a.s.l. (GH) to 29% at 450m a.s.l. 
(KL). Dominant species of the Zygoptera showed narrower trends compared to 
Anisoptera with most dominant species ranging from 18% at 790 m a.s.l. (BG) 
and 1050 m a.s.l. (CE) with Ceriagrion glabrum and Lestesplagiatus respectively 
to 23% at 100 m a. s. 1. (SB) and 1350 m a. s. 1. (GH) with Ischnura senegalensis 
and Pseudagrion kersteni , respectively. 

SITE VARIABLES AND ELEVATION 

ANOVA results (Tab. Ill) selected percentage shade (F = 53.52; P<0.001), float¬ 
ing and submerged vegetation (F = 20.43; P<0.001) marsh (F = 90.02; P<0.001), 
marginal forest (F = 335.39; P = 0.001), and marginal grasses (F = 15.56; P = 
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Table II 

Species of Odonata larvae recorded from all sites across the five elevations during the whole sampling period 
with species code names used in analyses 


Species 

Code 

Site abbreviations and elevation (m a.s.1.) 

SB KL BG CE GH 

(100 m) (450 m) (790 m) (1050 m) (1350 m) 

ZYGOPTERA 

Coenagrionidae 

Africallagma glaucum (Burmeister, 1839) 3 

Aglm 

L 

L 

L 

L 

L 

Agriocnemis falcifera Pinhey, 1959 4 

Afal 

L 

L 

L 

L 

L 

Ceriagrion glabrum (Burmeister, 1839) 3 

Cglm- 

L 

L 

L 

L 

- 

Ischnura senegalensis (Rambur, 1842) 3 

Isen 

L 

L 

L 

L 

L 

Pseudagrion kersteni (Gerstacker, 1869) 3 

Pker 

L 

L 

L 

L 

L 

P. massaicum Sjostedt, 1909 3 

Pmas 

L 

- 

L 

L 

- 

P. salisburyense , Ris, 1921 3 

Psal 

L 

L 

L 

L 

- 

Lestidae 

Lestes plagiatus (Burmeister, 1839) 3 

Lplg 

. 

L 

L 

L 

L 

L. tridens McLachlan, 1895 1 

Ltri 

L 

- 

- 

- 

- 

ANISOPTERA 

Aeshnidae 

Anax imperator Leach, 1815 3 

Aimp 

L 

L 

L 

L 

L 

Libellulidae 

Acisoma panorpoides Rambur, 1842 2 

Acis 

L 

. 

. 

L 

L 

Crocothemis erythraea (Brufie, 1832) 3 

Cery 

L 

L 

L 

L 

L 

Orthetrum caff rum (Burmeister, 1839) 3 

Ocaf- 

- 

- 

- 

L 

L 

O. Julia falsum Longfield, 1955 3 

Ojul 

L 

L 

L 

L 

- 

Trainee basilaris (Beauvois, 1817) 2 

Tbas 

L 

L 

- 

L 

L 

Trithemis dorsalis (Rambur, 1842) 3 

Tdor 

- 

L 

L 

L 

L 

T stictica (Burmeister, 1839) 1 

Tsti 

- 

L 

L 

L 

L 

T. arteriosa (Burmeister, 1839) 3 

Tart 

* 

- 

- 

L 

L 


1 Common African species whose range extends south just over the border into South Africa, but are 
local or rare in the country; 2 African species that are widespread and/or locally common in South 
Africa; 3 African species that are regularly seen in the right habitats; 4 Species endemic to South Af¬ 
rica i.e. south of the Limpopo River (Samways 2002). 

0.001) as showing highly significant differences across elevations. These results 
showed that habitats varied to a large extent across different elevations. 

LARVAL SPECIES RICHNESS, ABUNDANCE AND ELEVATION 

Larval Anisoptera species richness and abundance varied significantly positively 
with increasing elevation (SB) (larval species richness: F = 53.06; P = 0.000, larval 
abundance: F = 16.25; P = 0.004). Zygoptera species showed a generally negative 
correlation in richness with elevation, although this trend was not statistically sig- 
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nificant (larval species 
richness: F = 2.17; P = 

0.179), and a signifi¬ 
cantly negative corre¬ 
lation in abundance 
with elevation (larval 
abundance: F = 8.69; 

P = 0.049). 

SPEARMAN’S RANK 
CORRELATION OF 
VARIABLES AND 
LARVAE 

Anisoptera larval 
species richness was 
strongly positively 
correlated with eleva¬ 
tion (Elev: r = 0.894; 

P<0.01), with floating 
and submerged vege¬ 
tation (Fsveg: r = 0.9; 

P<0.001), with tur¬ 
bidity (Tur: r = 0.894; 

P<0.01), and with pH 
(r = 0.89; PO.01). An¬ 
isoptera larvae howev¬ 
er, were negatively cor¬ 
related with percentage shade (%Sh: r = -0.803; P<0.01). They showed a weak 
positive correlation with water temperature (Wt: r = 0.44; P<0.05) and marsh (r = 
0.57; P<0.05). Zygoptera larvae species richness was correlated with floating and 
submerged vegetation (Fsveg: r = 0.866; P<0.01), with marsh (r = 0.8; P<0.01), 
with turbidity (Tur: r = 0.87; P<0.01), and with pH (r = 0.8; P<0.01), although 
weakly correlated with water temperature (r = 0.5; P< 0.01). There was a weak 
negative correlation with percentage shade (%Sh: r = -0.44; P<0.01). These fac¬ 
tors appeared to be selected for at varying degrees by larval species (Tab. IV) i.e. 
elevation was more important for Anisoptera than for Zygoptera species. 

MULTIVARIATE ANALYSES 

Species-site-variable triplots for larvae of Anisoptera (Fig. 3) and Zygoptera 
(Fig. 4) showed that most assemblages related to certain site variables and indi- 



Fig. 2. Species accumulation curves for (a) Anisoptera, (b) Zygoptera 
larvae species using all 30 sub-sites. 
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Table III 

ANOVA results of measured site variables across all five elevations. Means followed by same letter 
across rows are not significantly different at 5% level of probability 


Site 

variables 

SB 

KL 

Sites 

BG 

CE 

GH 

F-ratio 

P-values 

%Sh 

15.25b a 

13.91b 

13.24b 

0.00a 

0.00a 

53.52 

<0.001 

At 

24.2 le 

26.26e 

26.14e 

27.06e 

27.06e 

1.85 

0.16 ns 

Fsv 

26.05d 

5.60c 

22.44d 

18.28d 

6.95c 

20.43 

<0.001 

Mfor 

30.0b 

47.5b 

25.12b 

0.00a 

0.00a 

335.39 

<0.001 

Mgra 

21.94 

7.39 

10.46 

19.50 

19.66 

15.56 

<0.001 

Marsh 

23.89g 

22.70g 

19.30g 

11.44f 

6.88f 

90.02 

<0.001 

pH 

7.013 

7.075 

7.088 

7.063 

7.041 

0.18 

0.948 ns 

Tur 

19.14 

25.76 

21.61 

19.86 

22.89 

1.60 

0.213 ns 

Wd 

22.30 

23.52 

23.31 

22.79 

21.21 

1.43 

0.26 ns 


SB= Stainbank Nature Reserve (100 m a.s.1), KL= Krantzkloof Nature Reserve (450 m a.s.1.), BG= 
Botanical Gardens (790 m a.s.1.), CE= Cedara (1050 m a.s.1.) and GH= Goodhope Estate (1350 m 
a.s.1.). 

cated how species responded or not to gradients of these variables. Important 
variables tended to be represented by longer arrows than less important ones. 
Accordingly, for the Anisoptera ordination, the first axis (horizontal) explained 
about 43.5% of the variation, with emergent and floating vegetation (EFveg), el¬ 
evation (Elev) and marginal grasses (Mgra) being most important. Less important 
variables were turbidity (Tur), and submerged vegetation (Sveg). The second axis 
(vertical) explained only 24.8% of the variation, with percentage shade (%Sh), 
marsh, water depth (Wd) and water temperature (Wt) being selected for. For the 
Zygoptera ordination, the first axis explained 40.6% of the total variation, with 
percentage shade (%Sh) being the main variable, and elevation less important in 
this axis. The second axis explained about 24.1%, of the total variation, select¬ 
ing for submerged vegetation (Sveg), marginal grasses (Mgra), water depth (Wd), 
marsh, turbidity (Tur) in decreasing order of importance. Both suborders showed 
great similarity with respect to broad biotope requirements, even though minor 
variations occurred with individual species preferences to site variables, which 
were usually synergistic with one another instead of having a singular effect. 

Anax imperator preferred emergent, floating and submerged vegetation condi¬ 
tions at biotopes. Most species preferred all conditions as explained by Anisop¬ 
tera species clustering at the centre of the ordination. Pseudagrion massaicum 
preferred increasing shady conditions while the remaining Zygoptera species also 
tolerated a combination of variables acting together hence are also clustered at 
the centre of the ordination. 

Intra-set correlations (Tab. V) and a summary of weightings attributed to the 
first two ordination axes for larvae of Anisoptera and Zygoptera (Tab. VI) indi¬ 
cated that species site variable correlations using CANOCO were strong. Cumu- 
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Table IV 

Spearman’s Rank Correlation coefficients between Anisoptera, Zygoptera larvae, odonate larvae rich¬ 
ness/abundance, and measured site variables across all elevations. *, 0.05>P>0.01; **, 0.01>P>0.001; 

***, P<0.001 


Site 

variables 

AniS 

AniN 

ZygS 

ZygN 

TOS 

TON 

Elev (m) 

0.894** 

0.900*** 

0.289 

-0110 

0.667** 

0.700** 

Pcir (m) 

0.112 

-0.600** 

-0.287 

-0.600** 

-0.154 

-0.20 

Mfor 

-0.229 

0.300 

0.148 

-0.335 

-0.026 

-0.205 

Fsveg 

0.900*** 

0.100 

0.866** 

0.800** 

0.564* 

0.600** 

Mgra 

0.205 

0.224 

0.00 

0.00 

0.205 

0.30 

Marsh 

0.574* 

0.100 

0.866*** 

-0.900** 

-0.667** 

-0.700** 

%Sh 

-0.803** 

-0.975*** 

-0.444* 

-0.051 

-0.658** 

-0.718** 

Wd 

-0.477* 

-0.800** 

-0.289 

0.200 

0.359 

-0.500* 

Tur 

0.894** 

0.500* 

0.866** 

0.700** 

.821*** 

0.900*** 

pH 

0.894** 

0.500* 

0.866** 

0.700** 

0.975*** 

0.900*** 

Wt 

0.447* 

0.900*** 

0.500* 

-0.500* 

0.447* 

-0.667** 


Elev = elevation, Pcir = reservoir circumference, Mfor = marginal forest, Fsveg = floating and sub¬ 
merged vegetation, Mgra = marginal grass, %Sh = percentage shade, Wd = water depth, Tur = tur¬ 
bidity, Wt = water temperature, Asp = AniS = Anisoptera species, AniN = Anisoptera individuals, 
ZygS = Zygoptera species, ZygN = Zygoptera individuals, TOS- total Odonata species, TON = to¬ 
tal Odonata individuals. 


lative percentage variance for larval species data and for larval species-environ¬ 
ment relations was 57.03% for Anisoptera and 59.4% for Zygoptera ordinations. 
This implies that measured site variables significantly accounted for the main 
variation in assemblage composition and distribution of larval species for both 
suborders. The respective eigenvalues, cumulative species variances and Monte 
Carlo tests for CCA are given in Table VI. 

DISCUSSION 

FAMILY AND SPECIES TRENDS 

The elevational tolerance of these species was remarkable. Six occurred at all 
elevations, 11 in two or more elevations, and only one species {Lestes tridens) oc¬ 
curring only at the low (SB) elevation. Even though these species were common, 
at least as adults, considerable sampling was required to detect them. Although 
species accumulation curves for both suborders were broadly similar, the rate at 
which new Zygoptera species appeared was higher than for Anisoptera, perhaps 
due to more microhabitat specificity in this suborder. 

The responses of pre-adult stages of temperate-centred species to low tem¬ 
peratures, show odonates to be principally tropical insects (NORLING, 1984), 
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Fig. 3. CCA triplot of Anisoptera Zygoptera larvae species (closed circles); 
sub-sites at elevations (open circles) and site variables (arrows) for pooled 
(2001 and 2002) data. Axis 1 is horizontal and axis 2 vertical. Species codes 
are as in Table II. Site abbreviations are: SB = Stainbank Nature Reserve 
(100 m), BG = Botanical Gardens (790 m), CE = Cedara (1050 m), and GH 
= Goodhope Estate (1350 m). Important site variables are: Elev = Eleva¬ 
tion; Mgra = marginal grasses; Sveg = submerged vegetation, Wd = water 
depth, %Sh = percentage shade. 


which decline in 
species richness 
with increasing 
elevation. In this 
study, Zygoptera 
species richness 
and abundance 
steadily declined 
with increasing el¬ 
evation from 100 
m to 1350 m, al¬ 
though not statis¬ 
tically significant¬ 
ly. Most species 
occurred through¬ 
out all elevations, 
with a few excep¬ 
tions e.g. Ceriagri- 
on glabrum which 
was not at 1350 
m and L. tridens 
which was only at 
100 m, being at its 
southern limit in 
this study. In con¬ 
trast to Zygoptera, 
Anisoptera larval 
species richness 


and abundance, as recorded here, increased significantly with increasing eleva¬ 
tion. This contrasts with the trend for adults where species richness declined with 
increasing elevation. Of the 47 adult species sampled during that study, larval 
stages of only 18 were identified here. Lack of suitable larval habitats, could have 
influenced larval species variations. While Trithemis stictica and Orthetrum caf- 
frum were easy to sample in the marshy shallows of pond edges at higher eleva¬ 
tions, this was not the case at low to middle elevations where the reservoir edges 
were steep, with little or no marshy shallow patches. Also, Africallagma spp. oc¬ 
curred in very shallow, marshy spots, could not be sampled as easily as open wa¬ 
ter. Anax speratus and Nesciothemis farinosa were rare as larvae and may have 
occurred among weed in the reservoirs. 
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SPECIES RICHNESS AND ENVIRONMENTAL VARIABLES 


Species richness in both odonate suborders had strong positive correlations with 
floating/submerged vegetation (Fsveg), pH and water temperature, while nega¬ 
tively correlated with shade conditions. Zygoptera species richness was positively 
correlated with marshy conditions e.g. Pseudagrion kersteni, P. salisburyense and 
R massaicum being found in these areas, whereas Anisoptera species richness 
was weakly correlated with this gradient. Larvae of Tramea basilaris and Anax 
imperator were mostly found in these marshy habitats. 

Species-site-variable triplots for both Anisoptera and Zygoptera indicated that 
as species clumped at the centres of the ordination diagrams, no particular meas¬ 
ured site variable clearly defined species assemblage distribution patterns. However, 
degrees of importance of variable gradients (denoted by length of arrows) could 
still be selected by using canonical ordination axes with their related intra-set cor¬ 
relation coefficients. Results showed that marginal grasses (Mgra), elevation (Elev) 


and emergent/sub¬ 
merged (EFveg) 
were the most 
important gradi¬ 
ents for Anisop¬ 
tera species, while 
percentage shade 
(%Sh) was very 
important for Zy¬ 
goptera. Other im¬ 
portant variables 
for both suborders 
were water depth 
(Wd), marsh and 
submerged vege¬ 
tation. 

Measured lo¬ 
cal variables e.g. 
vegetation (struc¬ 
tural and compo- 



- 1.0 + 1.0 


sitional), pH, wa¬ 
ter depth, turbid¬ 
ity and marshiness 
as well as regional 


Fig. 4. CCA triplot of Zygoptera larvae species (closed circles); sub-sites 
at elevations (open circles) and site variables (arrows) for pooled (2001 
and 2002) data. Axis 1 is horizontal and axis 2 vertical. Species codes are 
as in table 2. Site abbreviations are: SB = Stainbank Nature Reserve (100 


(large-scale) proc¬ 
esses e.g. climate 
(water and ambi- 


m), BG = Botanical Gardens (790 m), CE = Cedara (1050 m), and GH 
= Goodhope Estate (1350 m). Important site variables are: Elev = Eleva¬ 
tion; Mgra = marginal grasses; Sveg = submerged vegetation, Wd = water 
depth, %Sh = percentage shade. 
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Table V 

Intra-set correlations between each of the site variables and Canonical Correspondence Analy¬ 
sis axes for all Odonata larvae species (Anisoptera and Zygoptera) and site variables across all five 

elevations combined 


Anisoptera larvae Zygoptera larvae 

CCA axisl CCA axis 2 CCA axis 1 CCA axis 2 


Eigenvalues 

0.221 

0.126 

0.214 

0.128 

Sveg 

0.0812 

0.0867 

0.0304 

-0.5544 

EFveg 

0.6112 

0.0235 

-0.0561 

0.0016 

Mgra 

0.41050 

-0.3232 

0.2932 

-0.2045 

Marsh 

-0.2108 

-0.5634 

-0.0757 

-0.5174 

%Sh 

-0.4425 

0.4760 

0.4478 

0.3357 

Wd 

-0.3130 

-0.6026 

-0.1966 

-0.5729 

pH 

-0.0415 

-0.0707 

-0.1671 

-0.3081 

Wt 

0.1985 

0.1181 

0.2179 

0.1701 

Alt (Elev) m) 

0.5656 

-0.4348 

-0.4716 

-.0630 

Tur 

-0.3160 

-0.5543 

-0.1896 

-0.4388 


ent temperatures, elevation and reservoir size), explained variability in larval as¬ 
semblages (SANDIN, 2003; SANDIN & JOHNSON, 2002). In this study, meas¬ 
ured variables accounted for slightly over half of the variability in species assem¬ 
blage distribution patterns even though some of the variation could be accounted 
for by adult and larval behavioural patterns. 

ECOLOGICAL IMPORTANCE OF SITE VARIABLES TO LARVAE 

It is important for adult Odonata to select and colonise reservoirs where con¬ 
ditions for larvae are optimal, because larvae live longer and are often unable to 
move away from a habitat in times of adversity (OSBORN, 1995). WILDER- 
MUTH (1992) suggested that Odonata may react to visual cues when selecting for 
oviposition sites. Larger water areas are also of greater permanence than smaller 
areas, with the degree of permanence of a water body influencing the distribu¬ 
tion of larvae (CANNINGS & CANNINGS, 1987). 

Adult habitat preferences can be explained by a combination of thermoregu¬ 
latory requirements and selection of oviposition sites. Selection of micro-habi¬ 
tats suitable for thermoregulation is probably the most common form of ther¬ 
moregulation in ectothermic insects (MAY, 1976, 1984). The presence of larvae 
of Trithemis stictica , Crocothemis erythraea, Anax imperator, Africallagma glau - 
cum, Ischnura senegalensis and Ceriagrion glabrum in shallow marshes at reservoir 
edges at higher elevation (1050 m and 1350 m) probably showed that insolation 
strongly influenced growth and metamorphosis in some species. Larvae of Tra- 
mea basilaris and Pseudagrion massaicum, among other species, tolerated warm, 
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shady conditions at the 100 m site. Larvae occur where they obtain most of their 
preferred food, and be simultaneously protected from predators (CORBET, 1999). 
There were fish e.g. Tilapia sparrmanii in these reservoirs, that probably fed on 
the odonate larvae (SKELTON, 1993), especially in the shallows. There was per¬ 
haps also habitat and micro-habitat partitioning among larvae that may reflect 
past competition (OSBORN, 1995). 

Linked with food partitioning is substrate and vegetation, because prey would 
occur on specific substrates and vegetation. Most larvae were closely associated 
with submerged vegetation. This was observed in Trithemis spp., C. erythraea and 
some coenagrionids which clung onto partially submerged plants (OSBORN, 
1995). Submerged vegetation of the genus Chara , and Elodea densa, may have 
been selected to enable greater chance of finding prey. Also, although E densa 
is an invasive alien plant, it probably increased the complexity of local plant ar¬ 
chitecture (SAMWAYS, 1999). This, in turn, encourages establishment of cer¬ 
tain species. 


HABITAT PREFERENCE AND ELEVATION 

Elevation as a regional process alone was not as strong as its indirect influence 
on local environmental factors, especially on the Anisoptera larval assemblage. 
Zygoptera species distributions depended largely on local habitat, and indirectly 
on regional processes. However, odonate larvae in general, being the product of 


Table VI 

Weighting summaries of the first two axes of CA and CCA for larvae in terms of variances account¬ 
ed for by the two axes. Monte Carlo probability tests of significance are given for the first canonical 
axis (AX1) and all four axes. *P<0.05; **P<0.01 


Axes 

Weightings Anisoptera larvae Zygoptera larvae 

All 4 All 4 

CA CCA axes CA CCA axes 


AXES 

AX1 

AX2 

AX1 

AX2 

AX1 

AX2 

AX1 

AX2 

Eigenvalues 

0.264 

0.186 

0.221 

0.126 

0.314 

0.24 

0.214 

0.128 

SP-ENO 



0.932 

0.877 



0.835 

0.777 

CPVS 2 

32.2 

54.9 

27 

42.4 

32.4 

57.3 

22.2 

35.4 

CPVS-EN 3 



43.6 

68.4 



40.7 

65 

F-Ratio 



2.446 


5.557 


4.840 

2.033 

P-value 



0.005** 


0.005** 


0.03* 

0.005** 


1 Species-environment correlations; 2 Cumulative species variance of species data; 3 Cumulative spe¬ 
cies variance of species-environment relations 
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aerial dispersal in vagile adults, are restricted to particular habitats that they cannot 
change, even in adverse conditions. Since larvae sampled were mostly eurytopic 
species, distribution of species for both suborders therefore largely depended on 
corresponding dispersal of adults, most of which were elevation tolerant. Simi¬ 
larly, the larvae were also tolerant of a wide range of ecological conditions. 

Except for Agriocnemis falcifera , most species recorded here also occur across 
large areas of Africa, and many are common throughout South Africa (SAM¬ 
WAYS, 1999; 2002), although Lestes tridens is limited to warm coastal areas. The 
question then becomes: why are these species in this part of Africa so remark¬ 
ably tolerant of elevational changes, both as adults and larvae? The reason may 
lie with the fact that this is an El Nino Southern Oscillation (ENSO) area, which 
receives continual short-term and long-term climatic variation, with fluctuating 
water trends. This means that only the most habitat tolerant species (as larvae as 
well as adults) can survive in the long-term. 

CONCLUSION 

The larvae over this wide elevational gradient (1350 m) were mostly geographi¬ 
cally widespread habitat-tolerant species. There were a few exceptions (e.g. Lestes 
tridens, which is a low-elevation species). Even Agriocnemis falcifera (the only na¬ 
tional endemic species sampled here) was elevation-tolerant. Community compo¬ 
sition showed that there were no truly distinct larval assemblages at the various 
elevations. Despite the apparent diversity and heterogeneity of habitat types, the 
fine-scale micro-habitat conditions favoured by most larval species were remark¬ 
ably common to many of these habitats. While local environmental variables and 
behaviour were important in explaining variations in larval assemblage structure 
and taxon richness for the various reservoirs, there were no major differences 
between these reservoir sites at the regional level. The implication of this is that 
reservoirs in this part of Africa may be very rich in dragonfly species, but all of 
these species are geographically widespread and locally common. The upshot is 
that although these artificial water bodies are not increasing the ‘extent of oc¬ 
currence’, they play a major role in increasing ‘area of occupancy’. Whether this 
has genetic value for these already common species is uncertain. Furthermore, it 
appears that these species are remarkably habitat-tolerant and elevationally wide¬ 
spread through many millennia of impact and natural selection over many ENSO 
events. Survival is about vagility and colonization, with these artificial reservoirs 
being an ideal and relatively permanent local habitat for them. 

ACKNOWLEDGEMENTS 

The National Research Foundation and the Gay Langmuir Trust are thanked for financial sup¬ 
port. 



Elevation tolerance in African odonate larvae 


279 


REFERENCES 

BENKE, A.C., 1976. Dragonfly production an4 prey turnover. Ecology 57: 915-927. 

BENKE, A.C., 1978. Interactions among coexisting predators: a field experiment with dragonfly 
larvae. J Anim. Ecol 47: 335-350. 

BENKE, A.C. & S.S. BENKE, 1975. Comparative dynamics and life-histories of coexisting dragon¬ 
fly populations. Ecology 56: 302-317. 

BUCHWALD, R., 1992. Vegetation and dragonfly fauna — characteristics and examples of bioce- 
nological field studies. Vegetatio 101:99-107. 

CALDWELL, J.P., J.H. THORP & T.O. JERVEY, 1980. Predator-prey relationships among larval 
dragonflies, salamanders and frogs. Oecologia 46: 285-289. 

CANNINGS, R.A. & S.G. CANNINGS, 1987. The Odonata of some saline lakes in British Colum¬ 
bia, Canada: ecological distribution and zoogeography. Adv. Odonatol 3: 7-21. 

CARCHINI, G., M.J. SAMWAYS & P.M. CALDWELL, 1992. Description of ultimate instar lar¬ 
vae of five higher altitude Trithemis species in southern Africa (Anisoptera: Libellulidae). 
Odonatologica 21: 25-38. 

COLWELL, R.K. & J.A. CODDINGTON, 1994. Estimating terrestrial biodiversity through extra¬ 
polation. Phil Trans: R. Soc Lond. 345: 101-118. 

CORBET, P.S., 1957. The life history of the dragonfly, Anax imperator Leach (Odonata: Aeshni- 
dae). J. Anim. Ecol. 26: 1-69. 

CORBET, P.S., 1999. Dragonflies: behaviour and ecology of Odonata . Harley Books, Colchester. 

CROWLEY, P.H. & D.M. JOHNSON, 1982. Habitat and seasonality as niche axes in an odonate 
community. Ecology 62: 1064-1077. 

GOTELLI, N.J. & R.K. COLWELL, 2001. Quantifying biodiversity: procedures and pitfalls in the 
measurement and comparison of species richness. Ecol. Monogr. 4: 379-391. 

JOHNSON, D.M., R.N. BOHANAN, C.N. WATSON & T.H., MARTIN, 1984. Coexistence of 
Enallagma traviatum (Zygoptera: Coenagrionidae) in Bays Mountains Lake, Tennessee: an 
in situ enclosure experiment. Adv. Odonatol 2: 57-70. 

JOHNSON, D.M. & P.H. CROWLEY, 1980. Habitat and seasonal segregation among coexisting 
odonate larvae. Odonatologica 9: 273-308 

LAWTON J.H., 1970. A population study on larvae of the damselfly Pyrrhosoma nymphula (Sulzer) 
(Odonata: Zygoptera). Hydrobiologia 36: 33-52. 

MACAN, T.T., 1977. The fauna in the vegetation of a moorland fishpond as revealed by different 
methods of collecting. Hydrobiologia 55: 3-15 

MAY, M.L., 1976. Thermoregulation and adaptation to temperature in dragonflies (Odonata: Ani¬ 
soptera). Ecol Monogr. 46: 1-32. 

MAY, M.L., 1984. Energetics of adult Anisoptera, with special reference to feeding and reproductive 
behaviour. Adv. Odonatol. 2: 95-116. 

MORIN, P.J., 1984a. The impact of fish exclusion on the abundance and species composition of 
larval odonates: results of short-term experiments in a North Carolina farm pond. Ecol¬ 
ogy 65: 53-60. 

MORIN, P.J., 1984b. Odonate guild composition; experiments with colonization history and fish 
predation. Ecology 65: 1866-1873. 

McCAFFERTY, W.P., 1984. Aquatic entomology. The fishermans and ecologist’s illustrated guide to 
insect and their relatives. Jones & Bartlett, Boston. 

NORLING, U., 1984. Life history patterns in the northern expansion of dragonflies. Adv. Odona¬ 
tol 2: 127-156. 

OSBORN, R., 1995. Niche and life-history differences in five highly sympatric species of Trithemis drag¬ 
onflies (Odonata: Libellulidae). Unpublished MSc. thesis, Univ. Natal, Pietermaritzburg. 

PIERCE, C.L., P.H. CROWLEY & D.M. JOHNSON, 1985. Behaviour and ecological interactions 



280 


A.S. Niba & M.J. Samways 


of larval Odonata. Ecology 66: 1504-1512. 

SAMWAYS, M.J., 1999. Diversity and conservation status of South African dragonflies (Odonata). 
Odonatologica 28: 13-62. 

SAMWAYS, M.J., 2002. A strategy for national red listing invertebrates based on experiences with 
Odonata in South Africa. Afr. Ent. 10: 157-166. 

SAMWAYS, M.J., R. OSBORN & I. VAN HEERDEN, 1996. Distribution of benthic invertebrates at 
different depths in a shallow reservoir in the KwaZulu-Natal Midlands. Koedoe 39: 69-76. 

SANDIN, L., 2003. Benthic macroinvertebrates in Swedish streams: community structure, taxon 
richness, and environmental relations. Ecography 26: 269-282. 

SANDIN, L. & R.K. JOHNSON, 2002. Ecoregions and benthic macroinvertebrate assemblages of 
Swedish streams. JIN Am. benthol. Soc. 19:462-474. 

SKELTON, P., 1993. A complete guide to the freshwater fishes of southern Africa. Southern Book 
Publishers, South Africa. 

STERNBERG, K., 1994. Niche specialisation in dragonflies. Adv.Odonatol. 6:117-198. 

TER BRAAK, C.J.F., 1986. Canonical correspondence analysis: a new eigenvector technique for 
multivariate direct gradient analysis. Ecology 67:1167- 1179. 

TER BRAAK, C.J.F. & C.W.N. LOOMAN, 1995. Regression. In: R.H.G. Jongman, C.J.F. Ter Braak 
& O.F.R. Van Tongeren, [Eds], Data analysis in community and landscape ecology, pp. 29-77. 
Cambridge Univ. Press, Cambridge. 

THORP, J.H. & E. A. BERGEY, 1981. Field experiments on responses of a freshwater benthic micro- 
invertebrate community to vertebrate predators. Ecology 62: 365-375. 

THORP, J.H. & M. COTHRAN, 1984. Regulation of community structure at multiple intensities of 
dragonfly predation. Ecology 65:1546-1555. 

TRAVIS, J., W.H. KEEN & J. JULIANA, 1985. The role of relative body size in a predator-prey re¬ 
lationship between dragonfly naiads and larval anurans. Oikos 45: 59-65. 

UBUKATA, H., 1981. Survivorship curve and annual fluctuation in the size of emerging populations 
of Cordulia aenea amurensis Selys (Odonata: Corduliidae). Jap. J. Ecol. 31: 335-346. 

UGLAND, K.I., J.S. GRAY, K. ELLINGSEN & E. ELLINGSEN, 2003. Species-accumulation 
curve and estimation of species richness. J. Anim. Ecol. 72: 888-879. 

WARD, J.V., 1984. Ecological perspectives in the management of aquatic insect fauna. In: V.H. Resh 
& D.M. Rosenberg, [Eds], The ecology of aquatic insects , pp. 558-577. Praeger, New York. 

WISSINGER, S.A. 1988. Life history and size structure of larval dragonfly populations. Jl N. Am. 
benthol Soc. 7:13-28. 

WEIR, J.S., 1974. Odonata collected in and near seasonal pools in Wankie National Park, Rhode¬ 
sia, with notes on the physico-chemical environments in which nymphs were found. J. ent. 
Soc sthn Afr. 37: 135-145. 

WILDERMUTH, H., 1992. Visual and tactile stimuli in choice of oviposition substrates by the drag¬ 
onfly Perithemis mooma Kirby (Anisoptera: Libellulidae). Odonatologica 21: 309-321. 



